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Using potential models we analyze range corrections to the universal law dictated by the Efimov 
theory of three bosons. In the case of finite-range interactions we have observed that, at first order, 
it is necessary to supplement the theory with one finite-range parameter, U, for each specific n- 
level [Kievsky and Gattobigio, Phys. Rev. A 87, 052719 (2013)]. The value of r® depends on the 
way the potentials is changed to tune the scattering length toward the unitary limit. In this work 
we analyze a particular path in which the length rg = a — as, measuring the difference between 
the two-body scattering length a and the energy scattering length as, results almost constant. 

Analyzing systems with very different scales, as atomic or nuclear systems, we observe that the 
finite-range parameter remains almost constant along the path with a numerical value of Fq ~ 0.87 
for the ground state level. This observation suggests the possibility of constructing a single universal 
function that incorporate finite-range effects for this class of paths. The result is used to estimate 
the three-body parameter k * in the case of real atomic systems brought to the unitary limit thought 
a broad Feshbach resonances. Furthermore, we show that the finite-range parameter can be put in 
relation with the two-body contact C 2 at the unitary limit. 


I. INTRODUCTION 

The study of universal behavior in few-boson systems 
is an intense subject of research nowadays. Universal 
properties appear for example in loosely bound systems 
in which the particles stay most of the time outside the 
range of the interaction. In this situation the details 
of the interaction between components are not impor¬ 
tant and the dynamics of the system can be described in 
terms of few control parameters. In the two-body sys¬ 
tem the scattering length a can be used as the control 
parameter. In fact, when the two-body binding energy 
Ei —> 0 its value results E 2 « h 2 /ma 2 and all of the two- 
body observables can be written in terms a as well flj]. 
The paradigm for this universal behavior is represented 
by the zero-range theory: the particles are outside the 
interaction range all of the time. In this case the above 
relation becomes exact, E 2 = K 2 /ma 2 , and the usual ef¬ 
fective range expansion for the s-wave phase-shift reduces 
to fccott) = — 1/a (with the total energy E = K 2 k 2 /m). 

Application of this theory to the three-boson system 
produces the Thomas collapse [|j: the ground state en¬ 
ergy is unbounded from below. Besides this singular be¬ 
havior, the three-boson system shows a very peculiar be¬ 
havior as the two-body scattering length approaches the 
unitary limit, 1/a —> 0. In this limit, as has been shown 
by V. Efimov in a series of papers [3, 4|, a system of three 
identical bosons interacting through a two-body short- 
range interaction shows a geometrical series of bound 
states whose energies accumulate to zero. The ratio be¬ 
tween the energies of two consecutive states is constant 
and does not depend on the nature of the interaction. 
This particular behavior is known as the Efimov effect 
and its observation has triggered an enormous amount 


of experimental as well as theoretical work in different 
fields as molecular, atomic, nuclear and particle physics 
(for instance see Refs. 0 , 6 <] and references therein). 

The spectrum of the three-boson system close to the 
unitary limit is described by the Efimov equation (or Efi¬ 
mov radial law) which can be expressed in a parametric 
form as follow 


E%/(h 2 /ma 2 ) = tan 2 £ 


/c*a = e ( n - ra *W s o e 


-A(£)/2s 0 


cos£ 


(la) 
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with A(£) a universal function whose parametrization 
can be found in Ref. jlj and s 0 ~ 1.00624 is a universal 
number. The scale at which a particular set of eigenval¬ 
ues are selected among the infinite set of values is fixed 
by k*, called the three-body parameter, which defines the 
energy h 2 K 2 /m for n = n* at the unitary limit. Knowing 
the value of k* the spectrum in terms of a is completely 
determined. 

The ability of tuning a in atomic-trapped systems 
has allowed different experimental groups to measure 
the value of the scattering length a_ at which the 
three-body bound state disappears into the continuum 
(£ —>• — 7 r). From Eq. ([IJ we can see that measuring 
a_ = — e _A (’ I ’) //2s °/K* « — 1.50763/ft*, is an indirect way 
for measuring the three-body parameter «*. As an inter¬ 
esting result, it has been experimental found 0041 , and 
theoretically justified [Toj, EH’ that in the class of alkali 
atoms a-/t ss —9.5, with I the van der Waals length. 
Recently the same behavior has been seen in a gas of 
4 He atoms [ll . This fact has extended the discussion 
of universal behavior to analyzed the dependence of the 
three-body parameter on the two-body dynamics [ 13 ] . 



In the above discussion we have used that Eq. m and 
the parametrization of A(£) have been derived in the 
zero-range limit (scaling limit). On the other hand ex¬ 
periments and calculations made for real systems deal 
with finite-range interactions and, for this reason, finite- 
range corrections have to be considered Q- In Refs. (THI — 
0, the authors have solved the Schrodinger equation 
using potential models in order to observe in which man¬ 
ner finite-range corrections manifest in numerical calcu¬ 
lations. The following modifications to the Efimov radial 
law have been proposed 


a B 


E%/E 2 = tan 2 £ 
-A(£)/2s 0 

r 3 =-. 

n cos £ 


(2a) 
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In the above equation there are two different cor¬ 
rections: one comes from the two-body sector and is 
taken into account by substituting as for a, defined by 
E 2 = K 2 /ma 2 B , with E 2 the two-body binding energy if 
a > 0, or the two-body virtual-state energy in the oppo¬ 
site case, a < 0 [18] . The second correction enters as a 
shift T 3 in the control parameter k+clb- As we will show 
below the shift is almost constant close to the unitary 
limit. Moreover, the values of the three-body parameter 
and the shift depend on the energy level. 

In both Eqs. m and ED the three-body energy is ex¬ 
pressed as a function of the scattering length. In princi¬ 
ple, there can be different ways and different mechanisms 
to change the scattering length and tune it to a —>• ±oo, 
and we refer to a given protocol as a path to the unitary 
limit. In this work we explore one of these paths, or one 
class of paths, in which the length r B = a — as is kept 
constant. Along this path we obtain that the shift T 3 , 
which takes into account range corrections in the three- 
body sector, is independent, at a first order, on the kind 
of potential we use. This fact points to a kind of univer¬ 
sality of the first-order range corrections. 

One can speculate that this class of paths to the uni¬ 
tary limit is pertinent for the description of broad Fesh- 
bach resonances; using the universal corrections we can 
thus estimate the value of the 3-body parameter and 
compare it to what found in literature. Finally, we show 
that the shift Tq can be related to the two-body contact 
C 2 for a trimer at the unitary limit. 

The paper is organized as follow: in the Sec. [TT] we 
study the path to the unitary limit which describes a 
broad Feshbach resonance, and we analyze the behaviour 
of the two-body properties, in particular of rs• In Sec. Mil 
the three-body system is studied along that path using 
different potential models. We construct the Gaussian 
universal function, and we argue on its universality inside 
the potential models. In particular we show how to use it 
to predict the three-body parameter k* in real systems. 
In the same section we show the relation between the 
two-body contact and the shift. In Sec. IIVI we make our 
conclusions. 


II. PATH TO THE UNITARY LIMIT 


In the zero-range Efimov theory as = a, and this is the 
only way to move toward the unitary limit. On the con¬ 
trary, in a potential-based theory the finite-range char¬ 
acter of the interaction allows for different paths con¬ 
necting the physical point to the unitary limit. For ex¬ 
ample, in atomic-trapped systems, experimentalists use 
Feshbach resonance to modified the interatomic poten¬ 
tial. By changing the intensity of the magnetic field the 
two-body scattering length a is modified and for partic¬ 
ular values of the magnetic field it diverges. These reso¬ 
nances, which are further classified as broad or narrow, 
can be interpreted as a particular path along which the 
unitary limit is reached. Different theoretical descrip¬ 
tions of Feshbach resonances are available and the most 
common describes this process as a coupled channel sys¬ 
tem 0; still, broad resonances can be simply modeled 
by modifying the potential strength. In the following we 
analyze this option. We define 

V x (r) = XV (r ), (3) 

a potential with variable strength where V (r) is the po¬ 
tential that reproduces the binding energy E 2 and the 
two-body scattering length a of a particular system. The 
original potential corresponds to A = 1, and the unitary 
limit is reached decreasing the value of A down to a crit¬ 
ical value A c . Examples could be a dimer of two helium 
atoms, E 2 « 1.3 mK and a ~ 190 ao, or the deuteron, 
E 2 & 2.22 MeV and a ss 5.4 fm (here a is the triplet scat¬ 
tering length). To analyze these two cases we consider 
the LM2M2 potential of Aziz [20] for the helium dimer 
and, in the case of the deuteron, a combination of two 
Yukawians as the MTIII nucleon-nucleon potential [2lj . 
It should be noticed that these potentials reproduced a 
large set of two-body data and, for this reason, are con¬ 
sidered realistic potentials. On the other hand the Efi¬ 
mov radial law depends only on a or as, this suggests 
the use of this minimal information to construct a two- 
parameter potential as a Gaussian 

V (r) = U 0 e" T ' 2/r ° (4) 


with the strength and range fixed to describe the two 
experimental data. In the specific cases analyzed here 
the LM2M2 and MTIII interactions can be mimicked by 
a Gaussian of range ro = rj 1 ® = 10 ao and range r 0 = 
Tq N = 1.65 fm, respectively. 

Using the four different potential models, we have stud¬ 
ied the behavior of the energy-scattering length as, of the 
scattering length a, and of the effective range r e ff close 
to A c . In this region these two lengths are related by the 
effective range expansion as 


J_ _ 1 r eS 
as a 2 a 2 B 


( 5 ) 


If we define e = A — A c , we have a,a B = 0( 1/e), while 
r e fi = r u + 0(e), or r e s = r u + 0(1/a), if we want to 
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emphasize the dependency on the scattering length. For 
convenience we have introduced the length r u as the value 
of r e ff at the unitary limit. Formally r u can be obtained 
from the effective range expansion that, at that limit, 
takes the form 

lim C °^° = - Prlk 2 H- (6) 

a—too K Z 

where <5o is the s-wave phase-shift, k is the energy mo¬ 
mentum and P is the shape parameter. Therefore 

r u = lim r e ff = 2 lim 

a— J-oo k—>0 

What we have observed is that while r e s changes con¬ 
siderably from its value r u as the system moves away 
from the unitary limit, the length r B defined as 


lim 

a—too 


cotdn 


(7) 



r B = a - a B , (8) 

results almost constant, even if r B = r u j 2 + 0{ 1/a). To 
better analyze this fact we write Eq.([5j) in the form 

^ = i-T + 0 (i/“ 2 )- < 9 > 

If r B were strictly constant, we could replace in the above 
equation 2 r B by r u , its value at the unitary limit, result¬ 
ing in an universal relation for the ratio r e s/r u in terms 
of the variable r u /a 

— =1-0.5—. (10) 

r u a 

In Fig. |l]we show r e g/r u as a function of the inverse scat¬ 
tering length for different potential models. We can ob¬ 
serve that these very different potentials, describing phys¬ 
ical system with different scales as the atomic LM2M2 or 
the nuclear MTIII interactions collapse on a very narrow 
band; this is more evident for positive values of a where 
the physical points are located. The numerical results 
for the different values of A are given by the solid points 
whereas the solid lines are fits to the numerical results 
using a linear plus a quadratic r u /a term. The results for 
the two Gaussian potentials, given by the solid (green) 
circles lie on a single line. 

The behavior of r B is also shown in Fig. [I] in the case of 
the Gaussian interactions, by the (red) circles. The (red) 
dashed line is a fit to the numerical results using a linear 
plus a quadratic r u /a term. It should be noticed that in 
the expansion of 2r B /r u s=s l+Aor u / a +..., the coefficient 
of the linear term, Ao ~ —0.01 is more than one order 
of magnitude smaller than the corresponding one in the 
expansion of r e g/r u (which results to be always very close 
to 0.5) resulting in the almost constant behavior of r B . 
In the figure the position of the dimer formed by two 
helium atoms and the deuteron are explicitly shown. 


FIG. 1. (color online). The effective range r e g as a function of 
the inverse of the scattering length (in units of r„) for different 
potential models. The length 2 i'b in units of r u is also shown. 

III. THE THREE-BODY SECTOR 

The modifications to the Efimov radial law proposed in 
Eq.© have been derived using potential models follow¬ 
ing the path to the unitary limit discussed above. The 
potential strength has been varied in order to cover a 
wide range of values of a from positive to negative val¬ 
ues up to a_, the value at which the three-body system 
disappears into the three-body continuum. 

Let concentrate our analysis to the ground state of the 
three-body system, that means n = n* = 0, such that 
k* = Kq; in this case Eq. m and Eq. (l2bl) are related 
by 

a(£) ^ o-b(0 = L^/k* , (11) 

where we want to stress that a and a B are evaluated 
at the same angle £ and so they are not related to the 
same E 2 value but to the respective solutions of Eqs. (flail 
and (Eal) . The above relation extends the concept of al¬ 
most constant behavior of the difference a — a B from 
the two-body sector to the three-body sector. Defin¬ 
ing r* = Tq/k*, this length represents the distance be¬ 
tween th e zer o-range and the finite-range theory in the 
(1/a, — x/lEgl) plane and, if measured at fixed £ angles, 
it results (almost) constant. In order to further analyze 
this fact we discuss first the case of a Gaussian potential. 

A. Gaussian universal function 

0 G 

The three-boson ground-state energy, E 3 ’ , has been 
calculated with the two Gaussian potentials for different 
values of A. At A = A c the results verify 

K H e r He = K ^ N r gN o.488 (12) 
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where and k^ n are obtained from the energy at the 
unitary limit in both cases; thus, for a Gaussian of range 
Tq the ground-state binding energy at the unitary limit 
is Ef ss ?i 2 (r^/0.488) 2 /m. The ground state energies 
close to the unitary limit can be analyzed using Eq. d2]) . 
The results can be described with high accuracy using a 
shift Tq having a constant plus a l/(ft*as) term 


e 3 _ p3 

1 0 — 1 0,0 


r 3 
1 0,1 


(13) 


The parameters Tq q ~ 0.87 and Tg^ « —0.14 are in¬ 
dependent of the range of the Gaussian potential. The 
small value of Tg 1 guarantees an almost constant behav¬ 
ior of the shift along the path of variable strength. This 
behavior can be explicitly studied by casting Eq. m in 
the following form 


Tg e -A(«)/2s 0 

k*clb k»cib cos£ 


(14) 


The right hand side can be plotted against 1/k *clb look¬ 
ing for a linear behavior. This is shown in Fig. [2] where 
the solid points represent the calculations with the Gaus¬ 
sian potentials and the solid line is a fit using the expan¬ 
sion of Eq. 01. From the figure the almost constant 
behavior of Tg, close to the unitary limit, is confirmed. 
Moreover this analysis shows the completely equivalence 
between the different Gaussian potentials using k* as a 
scale factor. 



FIG. 2. (color online). Study of Ba. (ll4l) as a function of 
1 /k,8b. The solid circles are the numerical results for the 
r.h.s. of the equation. The (red) solid line is a fit to the 
calculations using a linear plus quadratic term. 


The equivalence between Gaussian potentials with dif¬ 
ferent ranges suggests the possibility of defining a Gaus¬ 
sian universal function A(£). The modified radial law of 


Eq.® can be written as 

Eg /E 2 = tan 2 £ (15a) 

e -A(£)/2so 

K.flfl = -T— , (15b) 

cos£ 

and A(£) can be extracted from the numerical solutions 
of the ground state energy, Eg, of the three-boson system 
calculated with a Gaussian interaction as 

^ = S ° ln ( h*Kl/m ) ' (16) 

The Gaussian universal function calculated in this way 
incorporates finite-range effects along the specific path 
used to reach the unitary limit. We would like to stress 
that exists only one Gaussian universal function indepen¬ 
dently of the range and strength of the Gaussian used to 
calculate it. In Fig. [3] the Gaussian function A(£) cal¬ 
culated using the two Gaussian potentials is given by 
the solid (green) circles. The two set of data completely 
overlap and cannot be distinguished from which Gaus¬ 
sian potential they have been calculated. The Gaussian 
universal function is compared to calculations using the 
LM2M2 interaction, (blue) squares, and the Yukawian 
MTIII potential, (black) triangles. For the sake of com¬ 
parison the zero-range universal function A(£) is shown 
in the figure by the solid (red) line. 

Interestingly, these interaction models, very different 
in scale and functional form, give rise to equivalent func¬ 
tions that collapse in a narrow band indicating that they 
have similar Tg 0 values. In fact for the LM2M2 potential 
the value Tg 0 « 0.82 is obtained whereas for the MTIII 
potential Tg 0 Ri 0.85. This means that, up to first order, 
the range corrections close to the unitary limit are al¬ 
most universal. Moreover, we can conclude that close to 
the unitary limit a two-parameter potential as a Gaussian 
contains the essential ingredients required by the dynam¬ 
ics. The different scales (here we have explored atomic 
and nuclear scales) are absorbed in the control parameter 
k*clb- The extension of the analysis to the first excited 
state, n = 1 level, produce a shift with a constant term 
of rf 0 « 0.08 for the potentials under consideration. 


B. Estimation of k* 

Although the equivalent functions are very similar for 
the Gaussian potential and for the realistic ones, differ¬ 
ences are evident when these potential are used to calcu¬ 
late energy levels in the tliree-boson system. The Gaus¬ 
sian adapted to reproduce the values of the LM2M2 in¬ 
teraction in the two-body system has to be supplemented 
with a three-body force in order to reproduce the corre¬ 
sponding e nerg ies in the three-boson system. For exam¬ 
ple, in Ref. [22| the low energy dynamics of three helium 
atoms has been described using a soft potential model 
consisting in a two-body Gaussian plus a Gaussian hy¬ 
percentral force adapted from the LM2M2 interaction. 
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FIG. 3. (color online). The zero-range universal function, 
(red) solid line, as a function of the angle £. The (green) 
circles, (blue) squares and (black) triangles are the equivalent 
functions calculated for the corresponding potentials. The 
dashed lines are the tangents at the unitary limit. 


Moving from the physical point towards the unitary limit 
this strategy implies that a set of calculations have to be 
performed in order to adapted the soft potential model 
to the realistic one. 

Instead of this type of approach here we want to use 
the Gaussian universal function to estimate the three- 
body parameter, «*, knowing a, £2 and £3 for some 
specific system. This can be done observing that Eq. m 
is a one-parameter equation as Eq. CD- The finite-range 
effects has been absorbed in A(£). An operative way 
to obtain this estimate is to determine the range and 
strength of the Gaussian from a and £ 2 . Calling the 
range determined in this way, the three-body parameter 
for that Gaussian can be determined from the relation 
of Ea. dT^l) . nf = 0.488/r^. Finally the three-body pa¬ 
rameter of the system can be obtained from the following 
relation derived using Eq. (fl5l) at a fixed value of the angle 
£ 


ifying £g /£ 2 = tan £ 2 = 97.0. Using the above equa¬ 
tion, we obtain an estimate of the three-body param¬ 
eter k* ft 0.0437aQ 1 to be compared to the value of 
0.0440aQ 1 obtained calculating the three-boson ground 
state energy with the LM2M2 potential using Eq. (J3]) with 
A c = 0.9743. In the case of the MTIII potential a sim¬ 
ilar reasoning produces and estimate of the three-body 
parameter of k* ft 0.269fm~ 1 compare to the exact cal¬ 
culation giving a value of 0.277 fm -1 . This results are 
well below a 3% accuracy, and a similar accuracy has 
been obtained for example using a next to the leading 
order effective field theory description [23| . 


C. Finite-range parameter in terms of the contact 


The shift modifies the Efimov radial law equation to 
take into account finite-range corrections [l~ 6 l j : thus, we 
have shown that the description of an Efimov-state en¬ 
ergy requires two three-body parameters, k* and the 
shift. We can relate the two parameters to other proper¬ 
ties of a bosonic gas at the unitary limit: the two- and 
three-body contacts [2l - (26| . 

We briefly recall that the two- and three-body contacts 
are defined by the following relations 



K 2 

87 Tma 


C 2 



m 


(18a) 

(18b) 


They are important quantities entering in several prop¬ 
erties of a many body system. For instance, the tail of 
the momentum distribution n(k) of a bosonic gas con¬ 
tains C 2 and C 3 in its asymptotic behaviour 


n ( k ) ^4 1 + ~ k ^° 3 ' ( 19 ) 

In Eq. (flTH) £(fc) is an universal log-periodic function [27|— 
l29j . it is not relevant for what follows. 

In Refs. [13, do] C '2 has been calculated for the Efimov 
trimer in the unitary limit, a —> ±00, and the result can 
be expressed in the following way 


K* 



(17) 


0 G 

where £ 3 ’ is three-boson ground state energy calculated 
using the Gaussian potential at a strength value such that 
E®/E 2 = £g’ G /£^ and with £^ the two-body energy 
calculated at that strength. 

For example, using the LM2M2 interaction the he¬ 
lium dimer and trimer energies are £2 = 1.303 mK and 
£3 = 126.4 mK respectively, and E 3 /E 2 = tan £ 2 = 97.0. 
Using the Gaussian representation with = 10ao 

and a strength Vo = 1.24294K the following values are 
obtained: £^ = 1.6209 mK and £g’ G = 157.3 mK, ver- 


C 2 = 8-7T — - k* ft 53.01 k* , (20) 

so 

where A'(— 7 r/ 2 ) ft 2.125850 (3C| is the derivative of the 
universal function calculated at the unitary point. In the 
same way, we can consider the two-body contact for the 
ground-state trimer calculated with the Gaussian poten¬ 
tial; using Eq. m we obtain 

C 2 = 8 tt A (~ 7T / 2 ' ) ^ 96.924 k* , (21) 

so 

where we have used the numerical derivative A' (-7T/2) ft 
3.8565. 
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In order to show how the two-body contacts, for the 
Efimov and the Gaussian trimers, are related to the shift 
Tq, we use Eqs. (l2bl) and (U5bll to write 

-A(f)/2s 0 -A(£)/2s 0 

rg = --— - --— ■ (22) 

COS 4 COS £ 

We can Taylor expand the right hand side of Eq. (E^l) 
around the unitary limit, £ = —n/2 obtaining 

= i - A '(-^/ 2 )) ■ (23) 

Therefore we conclude that the constant term of the shift 
is proportional to the difference of the derivatives of the 
Gaussian and zero-range universal function at the unitary 
limit. The derivatives are shown as the dashed lines in 
Fig. |3] Combining Eqs. (E51) . KU1) . and Oil) we obtain 

Fo o = 77 T -—(^2 — C 2 ) « 0. 86 5 3 . (24) 

’ lu7T ft* 

The shift can therefore be related to a different varia¬ 
tion of the three-boson energy with respect to the scat¬ 
tering length at the unitary limit given by the potential 
models with respect to the zero-range case. In the case 
of real potentials moving along the path in which tb is 
constant the shift measures the increase in this variation. 
Here we have observed almost equal values of the shift 
using different potential models with variable strengths 
to explore the dynamics close to the unitary limit. How¬ 
ever, we point out that in systems governed by a coupled- 
channel dynamics, different values of Tj) 0 are possible, 
even negative ones. 

In the case of the three-body contact, with the defini¬ 
tion given in Eq. fllSbj) it results C 3 = ft 2 in both cases, 
using the zero-range theory or the Gaussian model. As 
before we can related the latter to the range of the Gaus¬ 
sian as C 3 = 0.238/(r(f) 2 . 


D. Relation with experimental results. 

In atomic trapped systems the measured three-body 
parameter is a_ and not the binding energy at the uni¬ 
tary limit. Therefore the method described before used 
to predict the three-body parameter from a known en¬ 
ergy value cannot be applied. However, it is possible to 
estimate the three-body parameter ft* using the value 
of a-. In fact, Eq. m can be solved at £ = — 7 r re¬ 
ducing the strength of the Gaussian potential up to the 
point in which the three-body energy results to be zero. 
From detailed calculations around such a point we obtain 
A(— 7 r) « —1.83 and, accordingly, 

ft*a)) = —2.483 . (25) 

It should be noticed that the parametrization proposed 
in Eq. (fl3l) can be used in connection with Eq. (B^l) to ap¬ 
proximate A(— 7 r) using the value of the zero-range uni¬ 
versal function, A(— 7 r) = —0.8262, derived in Ref. [31,] 


and solving a second order equation. In such a case we 
obtain an estimate of ft*a^ ~ —2.44 in reasonable agree¬ 
ment with the above (more exact) result, showing that 
the two-ternr parametrization of the shift Tq can be ex¬ 
tended up to the three-body threshold into the contin¬ 
uum. Also the constant relation between the scattering 
length and the energy scattering length holds at that 
threshold, a^ = a_ — tb, and it can be used to deduce 
the three-body parameter in terms of a~. For systems 
having a van der Waals tail we propose the following es¬ 
timate 


-2.48 
a- - r B 


0.22 


(26) 


where we have introduced the van der Waals length £ 
and we have approximate the ratios a_/£ ~ —9.5 and 
7 \b/£ ~ 1.4. The first ratio reflects a class of universality 
determined by the van der Waals length (see Ref. [ill] 
and references therein) whereas the second ratio can be 
verified using theoretical models as has been done in the 
present work in the case of He or, in the case of Cs, by 
the model given in Ref. [II]. To be noticed that a very 
close relation to determine ft* in the case of two-body 
Leonnard-Jones potentials has been given in Ref. [351 ] 
showing that also this kind of potentials are well de¬ 
scribed using the Gaussian universal function. 

Though approximate, the estimate for the three-body 
parameter given in Eq. (l26l) agrees well with estimates 
given in the literature. For an helium trirner, using 
£ = 5.1 ao, Ea. fliTOl) predicts a three-body parameter of 
ft* ~ 0.043 a ,]" 1 whereas calculations using the LM2M2 
potential give ft* ~ 0.044ag 1 [13,EH- It is interesting to 
noticed that the above relation is also in agreement with 
the estimates of Ref. [34j for the three-body parameter 


of 6 Li (ft* = 0.00678 a^ 1 ) and 133 Cs (ft* = 0.0017ao). 


In the first case, using £ = 31.26 ao Eq. 
0.007a,” 1 

the second case, using £ = 101 ao a value of ft, 


predicts 


ft* ~ U.UUYa 0 in good agreement with the estimate. In 

0.002ao 1 

is obtained, in reasonable agreement with the quoted 
value in Ref. (HI. 


IV. CONCLUSIONS. 

We have analyzed the behavior of a three-boson sys¬ 
tem close to the unitary limit using potential models dif¬ 
fering in form and scale. Specifically we have used the 
LM2M2 helium-helium interaction, the MTIII interac¬ 
tion constructed to describe the s-wave nucleon-nucleon 
interaction in triplet state and two Gaussian representa¬ 
tions of these potentials. The strength of the potentials 
have been varied in order to move the system towards 
the unitary limit. A first observation was that along this 
path the length tb = a—as results to be almost constant. 
This quantity characterizes the system, it is about 7.2 a 0 
for the helium system and about 1.2 fm for the nuclear 
system. This can be understood as a particular path to 
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the unitary limit and the Gaussian representation of the 
potential is constructed to reproduce the length tb along 
that path. 

The analysis of the three-boson ground-state binding 
energy close to the unitary limit using the Gaussian po¬ 
tentials have shown different scaling properties allowing 
to determine the three-body parameter, from the 
range of the Gaussians. Moreover this family of poten¬ 
tials produces a shift with a constant term of Tq 0 ss 0.87, 
independently of the range, and a small first order cor¬ 
rection proportional to 1 /k*ob. This characteristic al¬ 
lows to construct a Gaussian universal function A(£). 
The shift has been shown to be related to the differ¬ 
ence between the zero-range-trimer two-body contact and 
the ground-state-trimer two-body contact at the unitary 
limit. Being the contact an observable, this connection 
allows for a measurement of the shift. Finally, the Gaus¬ 


sian universal function A(£) has been used to estimate 
the three-body parameter /c* knowing the two- and three- 
boson binding energies at a particular value of the two- 
body scattering length a. Since in atomic-trapped sys¬ 
tems estimates of the scattering length at dissociation are 
given, we make the analysis from the known value of er¬ 
as well. In this case we use the relations a-/l ~ —9.5 
and tb/1 ~ 1.4, verified in several atomic species. 

In many cases the zero-range universal function A(£) 
has been used to analyze experimental results with the 
conclusion that not always a complete agreement was 
found. In those cases in which the agreement is not sat¬ 
isfactory we suggest the use of the Gaussian universal 
function A(£) as an alternative to make the analysis. To 
this respect a parametrization of the function in terms of 
the angle £ could be helpful. Studies along this line are 
at present in progress as well as the analysis of different 
classes of paths to reach the unitary limit. 


[1] E. Braaten and H.-W. Hammer, 

Physics Reports 428, 259 (2006) 

[2] L. H. Thomas, Phys. Rev. 47, 903 (1935) 

[3] V. Efimov, Phys. Lett. B 33, 563 (1970) 

[4] V. Efimov, Sov. J. Nucl. Phys. 12, 589 (1971), [Yad. Fiz. 
12, 1080-1090 (1970)]. 

[5] F. Ferlaino, A. Zenesini, M. Berninger, B. Huang, H. C. 
Nagerl, and R. Grimm, Few-Body Syst. 51, 113 (2011) 

[6] T. Frederico, L. Tomio, A. Dolfino, M. R. Hadizadeh, 
and M. T. Yamashita, Few-Body Syst. 51, 87 (2011) 

[7] O. Machtey, Z. Shotan, N. Gross, and L. Khaykovich, 
Phys. Rev. Lett. 108, 210406 (2012) 

[8] S. Roy, M. Landini, A. Trenkwalder, G. Semeghini, 
G. Spagnolli, A. Simoni, M. Fattori, M. Inguscio, and 
G. Modugno, Phys. Rev. Lett. Ill, 053202 (2013) 

[9] P. Dyke, S. E. Pollack, and R. G. Hulet, 
Phys. Rev. A 88, 023625 (2013) 

[10] J. Wang, J. P. D’Incao, B. D. Esry, and C. H. Greene, 
Phys. Rev. Lett. 108, 263001 (2012) 

[11] P. Naidon, S. Endo, and M. Ueda, 

Phys. Rev. A 90, 022106 (2014) 


[12] S. Knoop, J. S. Borbely, W. Vassen, and S. J. J. M. F. 
Kokkelmans, Phys. Rev. A 86, 062705 (2012) 


[13] P. Naidon, E. Hiyama, 

and M. Ueda, 

Phys. Rev. A 86, 012502 (2012) 


[14] C. Ji, D. R. Phillips, 
EPL 92, 13003 (2010) 

[15] A. Kievsky and 

and L. Platter, 

M. Gattobigio, 

Phys. Rev. A 87, 052719 (2013) 


[16] M. Gattobigio and 

A. Kievsky, 

Phys. Rev. A 90, 012502 (2014) 



[17] A. Kievsky, N. K. Timofeyuk, and M. Gattobigio, 
Phys. Rev. A 90, 032504 (2014) 

[18] S. T. Ma, Rev. Mod. Phys. 25, 853 (1953) 

[19] C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, 
Rev. Mod. Phys. 82, 1225 (2010) 

[20] R. A. Aziz and M. J. Slaman, 
J. Chem. Phys. 94, 8047 (1991) 

[21] R. A. Malfliet and J. A. Tjon, 
Nuclear Physics A 127, 161 (1969) 

[22] A. Kievsky, E. Garrido, C. Romero-Redondo, and 

P. Barletta, Few-Body Syst. 51, 259 (2011) 

[23] E. Braaten and H.-W. Hammer, 

Phys. Rev. A 67, 042706 (2003) 

[24] S. Tan, Annals of Physics 323, 2952 (2008) 

[25] S. Tan, Annals of Physics 323, 2971 (2008) 

[26] S. Tan, Annals of Physics 323, 2987 (2008) 

[27] Y. Castin and F. Werner, 

Phys. Rev. A 83, 063614 (2011) 

[28] E. Braaten, D. Kang, and L. Platter, 

Phys. Rev. Lett. 106, 153005 (2011) 

[29] D. H. Smith, E. Braaten, D. Kang, and L. Platter, 
Phys. Rev. Lett. 112, 110402 (2014) 

[30] F. Werner and Y. Castin, (2010), arXiv: 1001.0774. 

[31] A. O. Gogolin, C. Mora, and R. Egger, 
Phys. Rev. Lett. 100, 140404 (2008) 

[32] M. Gattobigio, A. Kievsky, and M. Viviani, 
Phys. Rev. A 86, 042513 (2012) 

[33] E. Hiyama and M. Kamimura, 

Phys. Rev. A 90, 052514 (2014) 

[34] B. Huang, K. M. O’Hara, R. Grimm, J. M. Hutson, and 
D. S. Petrov, Phys. Rev. A 90, 043636 (2014) 

[35] D. Blume, Few-Body Syst. 56 , 859 (2015) 


7 





